Crossover from Fragile to Strong Glassy Behaviour in Kinetically Constrained Systems 
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We show the existence of fragile-to-strong transitions in kinetically constrained systems by study- 
ing the equilibrium and out-of-equilibrium dynamics of a generic constrained Ising spin chain which 
interpolates between the symmetric and fully asymmetric cases. We find that for large but finite 
asymmetry the model displays a crossover from fragile to strong glassy behaviour at finite tem- 
perature, which is controlled by the asymmetry parameter. The relaxation in the fragile region 
presents stretched exponential behaviour, with a temperature dependent stretching exponent which 
is predicted. Our results are confirmed by numerical simulations. 
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Glasses are everywhere in nature. It is difficult to find a 
liquid which when supercooled does not form the amor- 
phous, microscopically disordered solid we call a glass. 
The salient feature of supercooled liquids is their dramat- 
ical slowing down with decreasing temperature, signaled 
by the increase of their relaxation times and viscosities 
by several orders of magnitudes in a temperature range 
of a few decades. For reviews see 

Given the generic nature of the glassy state, universal 
principles for the classification of glass-forming materi- 
als are of paramount importance. Central to this is the 
concept of fragility |Q,[lJ , which measures the speed with 
which viscosity and relaxation times grow as a system ap- 
proaches the glass transition temperature. Liquids which 
display Arrhenius behaviour, that is, the logarithm of 
their viscosity or relaxation time grows linearly with in- 
verse temperature, are classified as 'strong' [§,0, as for 
example network liquids like Si02 and Gc02, and define 
one of the extremes of the classification by fragility. Most 
liquids, however, behave in a non- Arrhenius manner, and 
the larger the departure from Arrhenius behaviour, the 
more 'fragile' BJl], the most fragile liquids, which define 
the second extreme in fragility, being polymeric in nature. 
An exception to this classification is supercooled water 
H : at temperatures close to its melting point it behaves 
as extremely fragile, while near the glass transition it is 
very strong. Supercooled water has a fragile-to-strong 
transition. 

Some of the simplest systems which display the slow co- 
operative relaxation characteristic of glasses are the facil- 
itated kinetic Ising models, first introduced by Fredrick- 
son and Andersen jf|, in which glassiness is not a con- 
sequence of either disorder or frustration in the interac- 
tions, but of the presence of kinetic constraints in the 
dynamics of the system. Depending on whether the con- 
straints are isotropic || , or fully directed , these mod- 
els may behave as strong or fragile glasses, and the low 
temperature dynamics can be understood in terms of ac- 
tivation over energy barriers Jg-jljj]. They are particu- 
larly useful in the study of activated processes, which 
become highly relevant for supercooled liquids near the 
glass transition, but are not taken into account by ap- 



proximations like mode-coupling theory [ fl2[ or mean- 
field models HQ. 

The purpose of this Letter is to show that fragile-to- 
strong transitions also occur in kinetically constrained 
systems. We prove this for the simplest case of the one- 
spin facilitated Ising chain. We study a generalization of 
the kinetically constrained Ising chain which interpolates 
between the cases of symmetric and fully asymmetric 
constraints @. We show that for large but finite asym- 
metry the model displays a fragile-to-strong crossover at 
finite temperature. The crossover temperature is con- 
trolled by the asymmetry parameter, which also deter- 
mines the largest timescale and energy barrier of the 
problem. The relaxation in the fragile region presents 
stretched exponential behaviour, with a temperature de- 
pendent stretching exponent which is obtained analyti- 
cally. We performed extensive numerical simulations to 
confirm our results. 

We consider a chain of Ising spins a.- L 6 {0,1} 
(i = 1, . . . ,N), with periodic boundary conditions, and 
Hamiltonian H = J2i a i- The dynamics is restricted to 
single flips of spins which have at least one nearest neigh- 
bour in the up state. The rates for the possible transi- 
tions are in general different depending on whether the 
up neighbour is on the right or left, and are given by: 
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where b 6 [0,1] and e = exp(— 1/T). Detailed balance 
is obeyed, and the stationary distribution is the Boltz- 
mann distribution at temperature T for the Hamiltonian 
H. The parameter b sets the degree of asymmetry of 
the kinetic constraints. The limiting values b — 1/2 and 
6 = (or 1) correspond to the Fredrickson-Andersen 
(FA) model || and the asymmetrically constrained Ising 
chain (ACIC) Q], respectively. 

Due to the non-interacting nature of the Hamiltonian 
the thermodynamical properties of the model are triv- 
ial, and are the same for any b. The energy density is 
given by the concentration c of up spins (or 'defects'), 
which in equilibrium becomes c cq = e/(l + e). At low 
temperatures c is very small, and since defects facilitate 
the dynamics, the system slows down: isolated up spins 
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are locally stable and the system has to overcome energy 
barriers to evolve. 

In contrast with the statics, the low temperature dy- 
namics depends strongly on the value of 6 (except at T 
strictly zero fyl]). In the symmetric limit 6 — 1/2, which 
corresponds to the FA model, isolated defects can dif- 
fuse to the left (resp. right) by means of processes (ii) 
and (hi) [resp. (iv) and (i)] of Eq. ([!]). Each move re- 
quires the temporary creation of one defect, so there is 
a single activation barrier to diffusion AE = 1. This 
constancy of the energy barriers implies that relaxation 
times follow the Arrhenius law tfa ~ ex p(A-E/T), char- 
acteristic of strong glass behaviour PJll]|. The decay of 
the concentration of up spins, in the out of equilibrium 
regime, is well approximated by the coagulation process 
A + A — > A, which means that the typical lengthscale 
grows as I ~ (i/r) 1 / 2 ||. The situation in the asym- 
metric limit b = or 1, which corresponds to the ACIC 
model 0, is very different. Here the activated diffusion 
mechanism of the symmetric case is absent: a defect at 
a distance 2™" 1 < d < 2™ from the nearest defect (in the 
direction of the constraint) has to cross a barrier AE = n 
to move [ [To|] , i.e., barriers grow with the logarithm of the 
size of relaxing regions. This means that typical length- 
scales grow as I ~ t T ln 2 , which leads to a relaxation time 
at low temperatures of the form tacic ~ exp(l/T 2 ln2) 
Jlo[ . This is the Bassler law jfT5|,2 used as an alternative 
to the Vogel-Fulcher equation "J^[ to represent fragile 
behaviour. The absence of a finite temperature singu- 
larity is consistent with the trivial statics of the ACIC 
model. 

We now consider the behaviour at intermediate val- 
ues of the asymmetry b. We first focus on the relaxation 
towards equilibrium after a quench from infinite tempera- 
ture. When b is not far from the symmetric limit b = 1/2, 
the rates of Eq. (|l|) for reactions to the left and right are 
comparable and the symmetric diffusive mechanism is 
still effective: the behaviour is essentially that of the FA 
model. The region of large but finite asymmetry, that is, 
when b (or 1 — 6) is small,is more interesting. In this case 
rates (i) and (ii) of Eq. (|l]) are very much suppressed re- 
spect to rates (iii) and (iv) (or vice- versa), the system can 
only make use of the asymmetric mechanism to relax, and 
behaves like the ACIC. The timescales associated with 
the asymmetric process, however, grow as the relaxing 
regions become larger: the timescale for relaxation of a 



region of length 2 Tl 



< d < T 



is r 



(n) 



and the 



dynamics takes place in stages labeled by n [|10[. The 
timescale for the symmetric process, ~ [6(1 — 6)e]~ 1 , 
is large, but not infinite, except in the ACIC limit, and 
in contrast to the asymmetric one does not grow with in- 
creasing length. This means that eventually ts becomes 
comparable to for some value of n. This defines the 
last stage of asymmetric relaxation n* before the system 
switches to symmetric behaviour: 




n* ~ 1 -Tln[6(l - 6)]. 



(2) 
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FIG. 1. Defect concentration c as a function of time t 
after a quench from infinite temperature to T — 1/6, for val- 
ues of the asymmetry b = 1/2 (FA model), 10" 2 , 10" 3 , 10" 4 , 
10 -J , and (ACIC model). Inset: typical lengthscale I = 1/c 
against scaled time variable T ln t. The upper dashed line cor- 
responds to t 1 / 2 and the lower one to t Tln2 . Simulations were 
performed using a continuous time Monte Carlo/Metropolis 
algorithm [16, 17] for a system of 10 s spins. 



In Fig. [I] we show the decay of defect concentration c 
(the energy density of the system) as a function of time 
t after a quench from an initial state at Tq = oo to a low 
temperature T = 1/6, for several values of the asymme- 
try 6. For the 6 = ACIC limit, at least four plateaus 
in the concentration are visible, which correspond to the 
different stages in the dynamics. The number of plateaus 
decreases with increasing 6 in accordance with Eq. (0): 
there are three plateaus for 6 = 10~ 5 (n* = 2.9), two for 
6 = 1(T 4 and 1(T 3 (n* = 2.5 and 2.2), and only one for 
6 = 1(T 2 and 1/2 (n* = 1.7 and 1.1), the latter being 
the FA limit. The Inset shows how the typical length- 
scale I = 1/c grows with time. While in the ACIC case 
it follows i Tln2 (lower dashed line), for all the nonzero 
values of 6 the system eventually switches to the diffusive 
i 1 / 2 behaviour of the FA case (upper dashed line). 

Let us turn to the implications that a finite asym- 
metry has on the relaxation time of the system. For 
arbitrary 6, the symmetric and asymmetric processes 
compete. The relaxation timescale is given by t ~ 
(t^ 1 + t^jq) , assuming that the corresponding rates 
add up. The factor 6(1 — 6) in the symmetric relaxation 
timescale t$ can be interpreted as an entropic barrier: 
t s ~ exp{[A£ - Tln6(l - b)]/T}, where AE = 1. No- 
tice that this 'free energy' barrier is precisely n* of Eq. 
(g). Thus, the suppression of the symmetric mechanism 
due to 6 decreases with decreasing temperature. The con- 
sequence of this on the relaxation time for a fixed value 
of 6 is the following. At higher temperatures the asym- 
metric process dominates, r ~ tacic j and the system 
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displays fragile relaxation. At lower temperatures, the 
symmetric process becomes dominant, r ~ rg, and the 
behaviour is strong. The crossover temperature for this 
fragile-to-strong transition is determined by 6, 



1- v /l-41n[6(l-b)]/ln2 
2 ln[6(l - 6)] ' 



(3) 



corresponding to ts ~ tacic- 

In Fig. U we show the equilibrium relaxation time r 
as a function of inverse temperature 1/T. Given that 
the static properties of the system are known exactly 
it is simple to construct low temperature equilibrium 
configurations. This allows to study the equilibrium 
dynamics down to really low temperatures in contrast 
with most glassy systems where only out of equilibrium 
quantities are accessible. We obtain the relaxation time 
through the connected equilibrium autocorrelation func- 
tion, C(t) = N' 1 J2i( a i( t ) (T i( )) ~ c l v where r is de- 
fined by C(t) = e _1 C(0). Alternative definitions of the 
equilibrium relaxation timescale give similar results. The 
figure shows the following features: (i) in the 6 = 
ACIC limit, the relaxation time has the fragile behaviour 
tacic for all temperatures, as expected; (ii) for small 6, r 
crosses over from fragile tacic &t higher temperatures to 
strong ts behaviour at lower ones, the crossover taking 
place around T c given by Eq. (||) (e.g., 1/T C = 3.2 for 
b = 1Q- 5 and 1/T C = 2.9 for b = 1CT 4 ); (iii) for larger b 
the fragile region shrinks, and disappears completely at 
the FA limit b = 1/2. The displacement of the curves in 
the strong regime is due to the entropic barrier. Rescal- 
ing the relaxation time by r — » 6(1 — 6)r makes them 
collapse, as shown in the top-right panel of Fig. |[ The 
bottom-right panel gives the effective activation barrier 
AE = dlriT /d(l/T). The crossover here appears as a 
change from the linear growth in the ACIC to the con- 
stant barrier of the FA model. The crossover is sharper 
the smaller 6. The high temperature behaviour of r is 
exponential in 1/T, which gives the offset in the straight 
line for 6 = 0, and becomes irrelevant at low tempera- 
tures. The slope of AE is about 1.7 rather than 2/ In 2. 
This difference is due to the non-exponential nature of 
the autocorrelation, which does not affect the 1/T 2 be- 
haviour of the log of r. 

The fragile-to-strong crossover can also be observed 
in the behaviour of equilibrium dynamical quantities. 
The simplest one to study is the equilibrium persistence, 
Pit) = N~ 1 c~ 1 J2i(n t t>=o a i{t'))i which measures the 
fraction of defects of the initial configuration which have 
never flipped between times and t. It is closely related 
to the equilibrium autocorrelation C(i), but is free from 
the problem of the recurrence of defects which makes the 
analysis of the latter more tricky, particularly in one di- 
mension. 
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FIG. 2. Equilibrium relaxation time r as a function of in- 
verse temperature 1/T, for values of the asymmetry b — 1/2 
(FA model), 1CT 1 , 1(T 2 , 10" 3 , 1(T 4 , 1CT 5 , and (ACIC 
model). Top-right: rescaled time 6(1 — 6)r against 1/T. Bot- 
tom-right: effective activation barrier AE as a function of 
1/T. Simulations were performed for system sizes which var- 
ied from 10 5 for high temperatures to 5 x 10 6 for low ones, 
and data points were averaged over twenty runs. 



As mentioned before, in the asymmetric limit, the 
probability to flip a defect depends on the distance to the 
nearest defect in the direction of the constraint. In equi- 
librium, the probability distribution of these distances 
is independent of time. The persistence may be then 
approximated by the sum of the independent exponen- 
tial relaxation of defects at different distances from their 



neighbours, P(t) 



-t/r 



p n e 



-t/T 



where 



tq ~ 1 is the timescale associated with the initial T inde- 



\2" 



pendent transient, and p n = (1 — c eq ) — (1 — c cq )" is 
the probability for a defect to have a chain of spins of 
length 2 ,l_1 < d < 2™ next to it, which takes into account 
the fact that the relaxation time for the corresponding 

(n) 

distances is t a . The equilibrium condition is crucial to 
assume an independent relaxation of the different length- 
scales. If the initial configuration is an out of equilib- 
rium one, the probability distribution of distances evolves 
in time and the approximation considered is no longer 
valid. For short times the persistence is dominated by the 
fastest exponential decay, leading to — InP ~ t/f for low 
temperatures, where f = ro/c cq defines the timescale for 
short times. The long time behaviour may be estimated 
replacing the sum by an integral which is evaluated in 
the saddle point approximation, in a manner similar to 
that of Ref. |lq] . As a result we obtain a stretched ex- 

- (t/r AC icf ] 



ponential P(t) 
exponent 



exp 



= (1 + 1/T In 2)" 



with a stretching 



(4) 



Notice that this is an alternative method to obtain tacic 
to the one of Ref. jl(J . In the symmetric limit the persis- 
tence is simply the sum of the relaxation of defects with 
and without an up neighbour. At low temperatures it 
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reads P(t) ~ 2 c cq e 



-t/r 



(1 - 2c cq ) e-'/ Ts . The small 



time behaviour is similar to the asymmetric case, while 
the long time one is given by P(t) ~ exp(— t/rs). 

The behaviour of the persistence for b = 10 -5 and 
b = 10~ 3 at different temperatures is shown in Fig. |^. 
We present the data in a double log scale for P and a log 
scale for t to display the different stretching exponents. 
As in the case of the out of equilibrium relaxation, the 
decay is first dominated by the asymmetric process cor- 
responding to the fragile regime. For short times it is ex- 
ponential with a characteristic timescale f ~ e _1 , as de- 
scribed above. At longer times we see a change of slope in 
the plot, which corresponds to stretched exponential be- 
haviour, with an exponent given by Eq. (|4|). The stretch- 
ing region is larger for smaller b, as expected. The fragile- 
to-strong crossover then takes place, and the persistence 
becomes exponential again, now with a timescale t$. In 
the Insets we rescale time by a factor of e to superimpose 
the curves in the exponential regimes of short and long 
times. The two limiting lines correspond to exp(— t/f) 
and exp(— t/rs). 

We conclude with a comment on the explicit spatial 
asymmetry in the definition of the model studied here. 
It seems that to obtain other than strong behaviour in 
systems with kinetic constraints it is necessary to con- 
sider cases in which spatial isotropy is explicitly broken, 
like in the ACIC @ or its generalizations, which is a 
rather unphysical feature. This is also the case of sys- 
tems with interactions, but which display a dynamical 
behaviour similar to the spin facilitated models . The 
system considered in this work, however, can be defined 
in an alternative but explicitly spatially symmetric for- 
mulation, which also clarifies the relation between the 
asymmetry b and the timescale and lengthscale at which 
the fragile-to-strong crossover takes place [^0| . Consider 
& as a collective field b(t) which takes values and 1 
with a characteristic timescale for flipping t&. The sim- 
plest physical choice for this dynamics would be a Poisson 
process. b(t) may also have a spatial dependence, but on 
lengthscales larger than the typical ones for the spin sys- 
tem. Since (b(t)) = 1/2 spatial symmetry is unbroken. 
The instantaneous value of b(t) is either or 1, so for 
times smaller than 77, the behaviour is that of the asym- 
metric case. For times much larger than 77, the system 
effectively sees the average of b(t) and the behaviour is 
the symmetric one. 77, sets the timescale for the fragile- 
to-strong crossover. This argument is easily generalized 
to higher dimensions by considering a collective vector 
field instead. 
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FIG. 3. Persistence P as a function of time t, for b — 10 
(left panel) and b — 10 -3 (right panel). Notice that we plot 
— lnP(f) in a log-log scale. The dotted lines are the expected 
stretching exponents /3. Insets: same as main panels with 
time rescaled as t — * te. Details of simulations are the same 
as in Fig. Q. 
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